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V/-~j ' The Degasperis-Procesi equation with self-consistent sources(DPESCS) is derived. The Lax 

00 " 

representation of the DPESCS is presented. The conservation laws for DPESCS are constructed. 

The peakon solution of DPESCS is obtained by using the method of variation of constants. 
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1 Introduction 

Soliton equations with self-consistent sources (SESCS) have attracted much attention in recent 
years. They are important integrable models in many fields of physics, such as hydrodynamics, state 
physics, plasma physics, etc [T]-[TH] . For example, the KdV equation with self-consistent sources 
describes the interaction of long and short capillary-gravity waves [5 J . The KP equation with self- 
consistent sources describes the interaction of a long-wave with a short wave packet propagating 
on the x — y plane at some angle to each other [2] . The nonlinear Schrodinger equation with self- 
consistent sources represents the nonlinear interaction of an electrostatic high-frequency wave with 
the ion acoustic wave in a two component homogeneous plasma [6]. The constrained flows of soliton 
hierarchy may be regarded as the stationary systems of the corresponding integrable hierarchy with 
self-consistent sources |15|-|18|. Since the Lax representation for constrained flows can be deduced 
from the adjoint representation of the Lax representation of soliton equation p3] , there is a naturally 
way to find the zero-curvature representation for SESCS [15j-|18j. From this observation, the soliton 
equations with self-consistent sources may be viewed as integrable generalizations of the original 
soliton equations. A systematic way to construct the soliton equations with self-consistent sources 
is proposed in [T5]-|18|. 



The Camassa-Holm equation, which was implicitly contained in the class of multi-Hamiltonian 
systems introduced by Fuchssteiner and Fokas in [20J and explicitly derived as a shallow water wave 
equation by Camassa and Holm in [21], has the form 



Ut + 2wu x — u xx t + 3uu x = 2u x ii xx + uu xxx . (1-1) 

Since the works of Camassa and Holm, this equation has become a well-known example of integrable 
systems and has been studied from many different points of view. 



It is a natural question to ask whether there are other third order dispersive PDEs sharing 
the integrability properties of Camassa-Holm equation. An answer has been given in [19] where 
the method of asymptotic integrability was applied to a family of third-order dispersive PDE 
conservation laws 

Ut + c u x + ^u xxx - c?u xxt = (ciu 2 + c 2 u\ + C 3 UU XX ) X , (1.2) 

where a, cq, c\, C2, c% and 7 are some arbitrary constants. Only three equations in this family satisfy 
the asymptotic integrability conditions. They are the KdV equation, the Camassa-Holm equation 
and the following new equation 

ut + u x + 6uu x + u xxx - a 2 (u xxt + -u x u xx + -uu xxx ) = 0. (1.3) 

By a coordinate transformation, this new equation could be written as [22] 

Ut - u xxt + 4uu x = 3u x u xx + uu xxx . (1.4) 

This encourages us to study the family of equations [22] 

Ut - u xxt + {b + l)uu x = bu x u xx + uu xxx . (1.5) 

All the equations in this family have peakon solutions of the form 

N 

u = Y,Pj(-t) e ~ lX ~ qM , (1-6) 

i=i 

where pj and qj satisfy the dynamical system 



(u ->\d G N . dG N 

^ = - (6 - 1} ^-' * = -djT> (L7) 



the generating function Gn reading 

,\ 
2 



1 N 



j,k=i 



Let m = u — u xx . The equation (1.4) could be written as 

m t + m x u + 3mu x = 0, (1-9) 

which is called the Degasperis-Procesi (DP) equation. It is showed in [23] that both the Camassa- 
Holm and the DP equation are derived as members of a one-parameter family of asymptotic shallow 
water approximations to the Euler equations. They describe the unidirectional propagation of 
nonlinear shallow-water waves. The DP equation has a third order Lax pair and a bi-Hamiltonian 
structure. The existence of its global solutions are considered in [24J. A new integrable hierarchy 
was extended from the DP equation in [25] . In [26] , the N-peakon solutions of the DP equation have 
been derived. The N-soliton solutions of the DP equation are obtained in [27]. Interesting results 
on the solutions of this equation have also been obtained in |28J-|32j. 

The soliton equation with self-consistent sources was firstly studied by Melnikov in [1]- [3]. The 
problem of finding soliton solutions or other specific solutions for equations with self-consistent 
sources has been considered in the past by several authors [4]-[18j. The present paper falls in 
that line of research, aiming at deriving Degasperis-Procesi equation with self-consistent sources 
(DPESCS) and to find for its special explicit solutions. 

We first construct the high-order constrained flow of the DP equation. Based on it we estab- 
lish the DPESCS by regarding the constrained flow of DP equation as the stationary equation 
of DPESCS in the same way as in [15j-[18j. The Lax pair of the DPESCS is obtained, which 
means that the DPESCS is Lax integrable. In [33]- [34] we pointed out that soliton equations with 
self-consistent sources can be regarded as soliton equations with non-homogeneous terms, and ac- 
cordingly proposed to look for explicit solutions by using the methods of variation of constants. 
Applying this technique to DPESCS we have been able to find its peakon solutions and peakon- 



antipeakon solution. 

This paper is organized as follows. In section 2, we extend DP equation including self-consistent 
sources and construct its Lax representation. In section 3 we derive its conservation laws. In section 
4, the peakon and the peakon-antipeakon solution are obtained. In section 5, we mention some open 
problems. 

2 The DPESCS and its Lax pair 

2.1 The DPESCS 

First we construct the high-order constrained flows of the DP equation, then establish the DPESCS 
and describe how to derive the Lax representation for the DPESCS. 

It is known that the Lax pair for DP equation (1.9) is |22| 

i>xxx = ipx-mXip, (2.1a) 

1 2 

Ipt = -Jlpxx ~U1p x + (U X + oT)"0- ( 2 - lb ) 

Consider the following equations obtained from the spectral problem and its formal adjoint 
problem for n distinct real Xj. 

Qj,xxx = Qj,x ~ mXjqj, j = 1, • • • ,n, (2.2a) 

rj, X xx = r jtX + mXjVj, j = 1, • • ■ , n. (2.2b) 

It is not difficult to find that 

—L = -\. qj r h j = l,.-.,n. (2.3) 



It is known that the DP equation possesses a bi-hamiltonian structure [22], namely: 

m t = B 1 6 ^ = B 6 ^, (2.4) 

dm dm 



where 



B = m 2 / 3 d x m 1/3 {d x - dl)- l m l / z d x m 2 ^, (2.5) 

B 1 = d x {l-d 2 x ){A-d 2 x ), (2.6) 

H = —fmdx, (2.7) 

H-l = -- I u 3 dx. (2.8) 



The high-order constrained flow of the DP equation is obtained from (2.2) for n distinct Xj, 
requiring that the "potential" m obeys the following constraint 

Qj,xxx = Qj,x ~ mXjqj, (2.9b) 

Tj,xxx = Tj,x + rnXjTj j = 1, • ■ ■ , n, (2.9c) 

where atj, j = 1, . . . , n are arbitrary constants. 

According to the approach proposed in [15]- [IS], the DPESCS consists of the following equation 



— B ( SHl Y^ S ^ j ) 

5m ^ 3 5m 

and the equations (2.2), which by using (2.3) and taking ay = — | leads to the DPESCS 

1 - 
m t = —um x — 3u x m — - YJd(l — <9 2 )(4 — d 2 )(Xjqjrj), (2.10a) 

i=i 

%>« = qj,x-mXjqj, (2.10b) 

rj,xra = ^,x + mXjrj, j = 1, • • • ,ra. (2.10c) 



2.2 The Lax representation of the DPESCS 

Comparing the DPESCS to the DP equation, we may assume that the Lax presentation of the 
DPESCS has the form 

ipxxx = tpx-mXtp, (2.11a) 

1 2 

ipt = -ji'xx -uip x + (u x + —)ip + aip + bip x + ap xx , (2.11b) 

where a, b and c are some functions of q 3 and r 3 to be determined. Requiring that under (2.10b) and 
(2.10c) the compatibility condition of (2.11a) and (2.11b), namely ip X xxt = ifttxxxi leads to DPESCS 
(2.10a) enables us to find that 

o-x — O'xxx + 3b x m\ + bm x X + 3c xx m\ + 3c x m x X + cm xx X 

n 
= Xj2ocjd(l - 8 2 ){A - d 2 ){q 3 r 3 ), (2.12a) 

i=i 
— 3a xx — b xxx — 2b x — 3c xx + 3c x mX + 2cm x X = 0, (2.12b) 

3a x + 3b xx + 2c x + c xxx = 0. (2.12c) 

From (2.10b) and (2.10c) we obtain two identities 

3(1 - <9 2 )(4 - d 2 ){ qj rj) = -3Xj(rn x W(q 3 ,rj) + 3mW{q 3 ,r 3 ) x ), (2.13a) 

W(q 3 ,r 3 ) xxx - W(q 3 ,r 3 ) x = X 3 (2m x q 3 r 3 + 3m(q 3 r 3 ) x ), (2.13b) 

where W(q 3 ,r 3 ) = q 3 r 3tX — q 3 , x r 3 is the usual Wronskian determinant. 
From (2.12c) we obtain that 

3a + 36 x + 2c + Psx = 0. (2.14) 



Then (2.14) together with (2.12b) lead to 



1 3 

a-xxx ~ a x = -<9(1 - <9 2 )(4 - d 2 )c + (cm x X + -c x mX) x , (2.15a) 

2((6 + -c x )a; - (6 + -c x ) xxx ) = 2cm x X + Sc^mA. (2.15b) 

With the relations above we can rewrite (2.12a) as 

_I,9(1 - a 2 )(4 - 8 2 )c + X(3m(b + ^c^ + m x (b + ^)) 

n 

= A J>,d(l - d 2 )(4 - S 2 )^-). (2.16) 

i=i 

Here we suggest that 

c = ^A,-qf,T,-, (2.17a) 

i=i 

b+\c x = Y j B j W{q j ,r j ), (2.17b) 

i=i 

where Aj, Bj, j = 1, . . . ,n are some undetermined constants. 

Finally with some calculations to determine Aj and Bj we solve (2.12) for a, b, c yielding: 



J\l AA 2 

a = ) J 7 ,2 A 9 (3A(gjr J>:r - qj, xx rj) - kXjqjTj - 2Xj(qjrj) xx ), (2.18a) 

i=i •? 
n i AA 2 
b = Y^ ~2 A 2 -V ( X ( qjrj ' x ~ gj>r ^ + X i( q i r ^ x ^ (2.18b) 

" AA? 



In this way, we obtain the Lax pair for (2.10a) under(2.10b) and (2.10c) as follows 

?Pxxx = ipx - mXip, (2.19a) 

1 2 

El AA„- 
g A 2 _ A 2 (. 3X (lj r j,xx ~ Qj,x X rj) - 4A j <? J T j - 2X j (q j r j ) xx ))ip 
j=i i 

A 1 AA 2 
+CL ~2 A 2 -A 2 ^ 9j ' rj> ~ ^'^ + X i( q i r ^ x ^ x 

j=l 3 

j=l 3 

which means that the DPESCS is Lax integrable. 

3 The infinite set of conservation laws for the DPESCS 

With the help of the Lax representation of the DPESCS, we could find the conservation laws for 
the DPESCS by a well-known method. First we assume that m, u and its derivatives tend to 
when | a; | — > oo, and assume that qj, rj and its derivatives tends to when x — > — oo. Set 

r = % (3.i) 



then the identity 



d dhnp d dlnip^ 



dt dx dx dt 

together with (2.11) implies that DPESCS has the following conservation law: 

|(r) = J^) = §- x {u x + a-(u + b)T-{ l - + c){T x + T 2 )), (3.2) 

where a, b and c are given by (2.14). Here we define that Q = u x + a — (u + b)T — (t + c)(T x + T 2 ). 
Using (2.11a) gives rise to 

r - v xx = mX + 3rr x + r 3 . (3.3) 



We can have two kinds of expansions for T in the power series of A. The first expansion is in 
positive powers of A [22J, 

r = J>*A fc+1 , (3.4a) 

k=0 
oo 

n = X>A fc+1 . ( 3 - 4b ) 

k=0 

We note that the odd densities hik+i are exact derivatives. The first two densities are h$ = u, 
hi = u 3 , which yield the conserved quantities Ho and Hi in (2.7) and (2.8). The following densities 
are nonlocal because that would be an inverse of the operator (1 — o£) in the sequence. 

The second expansion could be 

oo 

r = ^/x fc A^, (3.5a) 

fc=0 

oo 

O = ^^A 1 ^. (3.5b) 

fc=0 

With the relations of (3.2) and (3.3), we can obtain the infinite densities and fluxes of the 
conservation laws. For brevity, we omit the recursion relations here. 

After some calculations, we can find that /ij. with odd subscripts are derivatives of some func- 
tions. So we define H^ s = J /jj2 S —2dx, which are taken as the conserved quantities. We could find 
the first few conserved quantities given by [Iq, \xi are as follows 

H_i = / m 1/3 dx, (3.6a) 



#_ 2 = — hm 2 x m^ 7/3 + 9m~ 1/3 )dx. (3.6b) 

The corresponding flux of the conservation laws are 

1 n 
G-i = m 1/3 (-u+ - ^2X 2 j (qjr j , x -q j , x r j )), (3.7a) 

3=1 

1 n 1 1 " 

G-2 = ~u x ^2 ^j(,Qj r j,xx ~ 9i,^ r i)) ~ -m x jn~ l {-u + - ^ A^r^ - qj, x Vj)). (3.7b) 



i=i 3=1 
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As the space part of the Lax pair is the same as that of DP equation, the densities of the 
conservation laws are also the same. Of course, as the time part is different, the fluxes will be also 
different. 

4 Solution of the DPESCS 

4.1 One peakon solution of the DPESCS 

As mentioned in the introduction, we will construct peakon solutions for DPESCS by the method 
of variation of constants. 

The DP equation has peakon solutions [22]. The one peakon is 

u = ce-' x - ct+a l, (4.1) 

where a is an arbitrary constant. The corresponding eigenfunction of (4.1) is 

q = p[sgn(x - ct + a){e~\ x - ct+a \ - 1) - 1], (4.2a) 

r = [3[sgn{x - ct + a){e-\ x - ct+a \ -1) + 1], (4.2b) 

where (5 is arbitrary constant as well. 

Taking a and (5 in (4.1) and (4.2) to be time-dependent a(t) and (3{t) and requiring that 

u = ce-lx-ct+aft)^ (43a) 

q = /3{t)[sgn{x - ct + a(t))( e - |a '- ct+Q(t)l - 1) - 1], (4.3b) 

r = p{t)[sgn{x - ct + a(t))( e - |x - ct+Q(t)l - 1) + 1] (4.3c) 

satisfies the DPESCS (2.10) for n = 1. We find that c = j-, a(t) can be an arbitrary function of t 

11 



and (3{t) = y/a!{t)c. So we have the one peakon solution for (2.9) with n = 1, Ai = A = \ 

u = ce -l*- c *+«MI, (4.4a) 

q = c ^c7(t)[sgn(x - ct + a(t))(e- lx ~ ct+a(t)l - 1) - 1], (4.4b) 

r = cy^At)[sgn(x-ct + a(t))(e- lx - ct+a{t)l -l) + l}. (4.4c) 

The one peakon of the DPESCS also has a cusp at its peak, located at x = ct — a(t). We note 
that, while for the one peakon solution of the DP equation travels with speed c and has a cusp at 
its peak of height c, for the DPESCS, the cusp is still at its peak of height c, but the speed of the 
wave is no longer a constant. 

4.2 The peakon-antipeakon solution of the DPESCS 

The DP equation has peakon-antipeakon solutions [22] 



u = SW*)C (e -|, +c |t|-H _ g-lx-cM+aK (45) 

where a is an arbitrary constant. The corresponding eigenfunctions of (4.5) is 



f3j -^—^ m - ) ( S9 n(x + c\t\ + a)(e-l*+ c l*l+"l - 1) 

-sgn(x - c\t\ + a)(e"l :r - c l t l +a l - 1)), (4.6a) 



c 



PJ {l _^_ 2cltl) (sgn(x + c\t\ + a)(e-l*+ c l*l+«l - 1) 

+sgn(x - c\t\ + a)(e" |:r - c|t|+a| - 1)). (4.6b) 



Taking a and f3 in (4.5) and (4.6) to be time-dependent a(t) and j3{t), and with the method of 
the variation of constants, the peakon-antipeakon solution of the DPESCS with n = 1 and Ai = ~ 



12 



IS 



U = ^_^ e -|z+c|t|+a W | _ e -|x-c|t|+a W |^ uj, 

1 — e — I I 



9 = C \! ~°i-%^ (°9n{x + c\t\ + a(t))( e -l^«l+«(«)l - 1) 

-s<7n(x - c|t| + a(*))(e- |a; - c|t|+a(t)l - 1)), (4.7b) 



r = cJ^^(^(x + c|t|+a(t))(e-l-H*l+«WI-l) 

+s<?n(x - c|t| + a(t))(e-l x - c l*l +a WI - 1)), (4.7c) 

where a(t) is an arbitrary function of t. 

5 Conclusion 

As a concluding remark, we want to stress that the construction of N peakon solutions for the 
DPESCS is still an open problem, as well as for the case of the Camassa-Holm equation with self- 
consistent sources. In fact, in the case of the Camassa-Holm equation and DP equation a reciprocal 
transformation is used, but so far we have not been able to extend it in the case of those equations 
with self-consistent sources. 
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